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ABSTRACT 

The albedo of snow for different cloudiness conditions is an important parameter in the Earth’s 
radiation budget analysis and in the study of snowpack thermal conditions. An efficient method is 
presented for approximate calculation of incident spectral solar flux and snowcover albedo in terms 
of atmospheric, cloud, and snow parameters. The global flux under partially cloudy skies is expressed 
in terms of the clear sky flux and a coefficient which models the effect of scattering and absorption 
by cloud patches and multiple reflections between the cloud base and snowcover. The direct and 
diffuse components of the clear sky flux are obtained using the spectral flux outside the atmosphere 
and the spectral transmission coefficients for absorption and scattering by molecules and aerosols. 

The high reflectivity of snow is considered when calculating the diffuse component. An approximate 
solution of the radiative transfer equation and parameterized water vapor attenuation coefficients are 
used for scattering and absorption by clouds. The spectral snow reflectance model consideres specular 
surface reflection and volumetric multiple scattering. The surface reflection is calculated using a crys- 
tal shape dependent bidirectional reflectance distribution function; volumetric multiple scattering is 
calculated using crystal size dependent approximate solution of the radiative transfer equation. The 
input parameters of the model are atmospheric precipitable water, ozone content, turbidity, cloud 
optical thickness, size and shape of ice crystals of snow and surface pressure. The model outputs 
spectral and integrated solar flux and snow reflectance as a function of solar elevation and fractional 
cloudcover. 
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The model is illustrated using representative parameters for the Antarctic coastal regions. The 
albedo for a clear sky depends inversely on the solar elevation. At high elevations the albedo depends 
primarily upon the grain size; at low elevation this dependence is on grain size and shape. The gradient 
of the albedo-elevation curve increases as the grains get larger and faceted. The albedo for a dense 
overcaste is few percent higher than the clear sky albedo at high elevations. A simple relation between 
the grain size and the overcaste albedo is obtained. For a set of grain size and shape, the albedo 
matrices (the albedo as a function of solar elevation and fractional cloudcover) are tabulated. 
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THE ALBEDO OF SNOW FOR PARTIALLY CLOUDY SKIES 


1. INTRODUCTION 

The albedo of snow is an important parameter in hydrology (Range, 1 975 and Coibeck cl al. 
1979) and climatology (Kukla and Kukla, 1974; Williams, 1975; Kellogg, 1975 and Ham. a;'.; Shukla, 
1976). The albedo of snow determines the total short wave radiation absorbed and refleued by snow. 
Die melting of snow as a result of solar heating is studied in hydrology for water resources manage- 
ment. Runoff from snow provides greater than 65% of the streamflow in the Western United States. 

The short wave radiation reflected by snow is needed for the Earth’s radiation balance analysis 
(Curran et. al. 1978). Predictions of global circulation models are sensitive to the extent of snow- 
cover because of its high reflectivity. A model for the albedo of snow can supplement airborne and 
spacecraft remote sensing observations using visible and near-infrared radiometers. 

The albedo of snow is defined as the ratio of reflected to incident solar energy. Calculation of 
the reflected energy requires using models for the incident spectral flux and the spectral reflectance 
of snow. Incident spectral flux varies with solar elevation, atmospheric parameters (e.g., precipitable 
water and turbidity etc.), the cloud characteristics (e.g., thickness and distribution of patches) and the 
reflectivity of snow. The spectral reflectivity of snow depends upon solar elevation, the grain size and 
grain shape, and the thickness of the snowcover. 

An approximate model for calculating the incident spectral flux for partially cloudy skies is dis- 
cussed in Section 2. The spectral integrated values (the insolation) are obtained and compared with 
observations. Section 3 discusses the spectral reflectance model of snow. In Section 4 the albedo 
values are obtained and compared with observations. 

2. PARAMETERIZATION OF INCIDENT SPECTRAL SOLAR FLUX 

A partially cloudy sky which includes the limiting cases of a clear or a totally overcast sky, is an 
inevitable reality of our weather systems. Clouds which are aerosols with diverse water phases, affect 
the spectral and directional characteristics of incident solar radiation. Clouds of different physical 
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characteristics (e.g. density and thickness), however, affect the incident radiation differently. The 
incident radiation for a totally cloudy sky is partially diffuse for high altitude thin clouds but is gen- 
erally totally diffuse for low altitude thick clouds. Also, near-infrared radiation is generally attenuated 
due to water vapor absorption. In a partially cloudy sky, the distribution of cloud patches is another 
factor which affects the incident radiation. 

In a clear sky the incident radiation consists of direct and atmospherically scattered (diffuse) 
components. Snow being highly reflective, a significant portion of this radiation is reflected back 
into space. Clouds over snowfields entrap this radiation, causing multiple reflection between the 
cloudbase and the snowcover. Because of spectral variation in the reflectivity of snow, this multiple 
reflection modifies the spectral characteristics of the incident flux. 

Accurate ' alculation of the incident flux is difficult for totally cloudy skies (Dave and Braslau, 
1975) and it is significantly more difficult for a partially cloudy sky (McKee and Cox, 1974). In 
many hydrologic and climatic studies, such accuracy is generally compromised for the sake of sim- 
plicity and efficiency by using parameterized equations of oiMerved or accurately calculated insola- 
tions in terms of fractional eloudeover (Sivkov, 1971; Lacis and Hansen, 1974). In the spirit of these 
parameterization equations, this section outlines a model for approximate calculations of the spectral 
fluxes, from which the insolation values can be obtained. The insolation results are compared with 
observations. Note that an understanding of the optical properties of snow (reflection, transmission 
and radiative heating) should be based upon the spectral values. 

2-1 . Parameterization Equations 

Fora partially cloudy sky. the incident flux may be written as (Wiscombe, 1975; Schneider and 

G o (a) [1 - c + cr, t (a)] 

C. An) -- — 1 — — (2-1 ) 


Dickinson, 1 976): 


where 


G 0 (a) * clear sky global flux for solar elevation a 
c * fractional cloudcover 
r eS (a) - transmission coefficient of cloud 
A cg * diffuse reflectance of cloudbase 
A f * diffuse reflectance of snow 

r A ■ transmission coefficient between cloudbase and snowcover. 

Noting that the clear sky fraction is (1 - c), the numerator of Equation 2*1 represents the sum 
of fluxes from both clear and cloudy portions of the sky. This numerator is essentially the basis of 
the Angstrom-Savinov parameterization of observed insolation in terms of fractional cloudcover 
(Sivkov, 1971). The denominator of this equation represents the enhancement of intensity due to 
multiple reflection between cloudbase and the snowcover (see Liljequist, 1956; Vowinckel and 
Orvig, 1962). 

For a partially cloudy sky, the instantaneous values of the flux are quite variable due to the 
dynamical nature of the position of solar disk with respect to the cloud patches. Therefore in data 
analysis and parameterizations time averaged (hourly, daily, or monthly) flux values are considered 
(Sivkov, 1971). The fractional cloudcover and the solar elevation are the average values correspond- 
ing to this time duration. The clear sky flux (G 0 ) is the sum of direct and atmospherically scattered 
components. The direct flux (Fj ) is calculated by multiplying the spectral flux outside the atmos- 
phere (F 0 ) by the parameterized spectral transmission coefficients for absorption by ozone (r 0 ), 
water vapor,(r H ), mixed gases (r f ) and scattering by molecules (r R ) and aerosols (r ( ), as outlined 
by Leckner (1978), from the equation 

F i = F o r o 3 t h 2 o r t r R T . si™ (2-2) 

where a is the solar elevation angle. The diffuse flux incident over snowcover arises from atmospheric 
scattering of direct flux and atmo jpheric backscattering of radiation reflected by snow. Because of 
the high reflectivity of snow, the diffuse flux incident over snowcover is generally larger than over 
other surfaces for similar atmospheric conditions. The diffuse flux (F a ) is calculated by the method 
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suggested by Robinson (1966) from the equation 


f 2 = * F o T o 3 n. 2 oV 


sin 4/3 a 


where the coefficient £ depends primarily upon the albedo of snow. For a clean fine-grained snow of 
high latitudes with albedo generally larger than 0.80 (c.g., Antarctica) this coefficient is 

£ = 0.75 for a < 10° 

0.70 + 0.005a for 10° < a < 40° 

The sum of fluxes (2-2) and (2-3) gives the clear sky global flux G 0 . All transmission coefficients in 
these fluxes are tabulated by Leckner (1 978). The cloud transmission coefficient (r cg ) can be written 
due to scattering by cloud drops (T $ ), and absorption by ambient water vapor (T A ) (Feigelson, 1978). 
For visible and near-infrared radiations, the cloud drops are almost pure scatterer of radiation with 
strongly asymmetric phase function. The delta-Fddington approximation of the radiative transfer 
equation (Joseph et al. 1976), which is designed for strongly asymmetric scattering medium, gives 
the scattering transmission coefficient as 


T* = e 


(2 + 3p o ) (1 - e” T *^°) -3r* (1 -g*)e' T /m ° 
4 + 3 (1 - g*) t * 


where 


t * = ( 1 - g 2 ) r 
g* = g/( 1 + g) 

= sina 

t = optical thickness of cloud 

g = scattering asymmetry factor (=0.8; Feigelson, 1978) 

Apart from scattering, a spectrally selective absorption occurs due to ambient water vapor. It 
is this absorption which produces meteorologically important radiative heating of clouds. In this 
study the following angular averaged transmission coefficient denved from Leckner ( 1978) is used: 


T* = cxp - 


0.498 X k 


(1 +41.915 X k) 0 ' 45 
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where 

X * water vapor path length within cloud 
k * spectral absorption coefficient of water vapor (Leckner, 1978). 

Although cloud drops are not pure scatterers (Zudnkowski & Strand, 1969; Herman, 1977), absorp- 
tion by Uquid water takes place roughly in the same spectral interval as for the gas. It is thus a good 
approximation to include the effect of absorption by cloud drops by adjusting the water vapor path 
length (Lads and Hansen, 1974). A mean water vapor path length for stratiform clouds is about 0.3 
cm (Lads and Hansen, 1974; Wiscombe, 1975). 

During the downward pass through the cloud, the strongly absorbing frequencies of the radia- 
tion are depleted by the cloud. As a result, the cloud is relatively transparent for the radiation which 
will be reflected by snow. For calculating the albedo of cloud (A ct ) for this upwelling diffuse radia- 
tion, it is a good approximation to treat the cloud as a pure scattering medium (Lads and Hansen, 
1974). The diffuse reflectance of cloud can be obtained from the scattering transmission coefficient 
(2-4) as: 

A c«'»-2 (Vo dM 0 

Jo 

2 ^1 e-'/?'* +3(1 -g*)T* 

% ■ ■ " * ■' ■ ■ ^ i ■ (2-6) 

4+3(1 -g*)r* 

For transmission coefficient between cloudbase and snowcover (r A ), one should consider aero- 
sols and water vapor attenuations because of their low scale heights. Assuming aerosols to be pure 
scatterers, its effect will primarily be to diffuse the radiation during the multiple reflection process. 
Since clouds are effective diffusor of the radiation, further scattering by the aerosol may be neglected 
in approximate calculations. Therefore, in this model only the selective spectral attenuation by water 
vapor within the cloudbase and snowcover is considered, i.e., the transmission coefficient (2-5) is sub- 
stituted for r A • The water vapor path length is taken to be one fourth the precipitable water in the 
clear atmosphere, which is typical for stratiform clouds (Wiscombe, 1975; Schmetz and Raschke. 
1979). 
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Note that in considering the effect of clouds on the incident radiation, the most important 
parameter is cloud's optical thickness because the reflection and transmission coefficients depend 


strongly upon it. Comparatively, the effect of moderate variations from the mean values of the water 
vapor path length within the cloud, and between cloudbase and snowcover is weak. 

Hie reflectance of snow shows strong spectral dependence, and the magnitude at different wave- 
length vary with the metamorphic changes (fcjellor, 1977). A recent study on the diffuse reflectance 
of snow based on the two-stream approximation of the radiative transfer equation (Choudhury and 
Chang, 1979) which provided qualitative understanding of these spectral and metamorphic variations, 
gave the following equation for a deep homogeneous snowcover: 

2(1 -gj) 1/2 


A* ” 1“ 


(1 - coa) ,/2 + (1 - to) 1/2 


(2-7) 


where 


w = ^ exp (-1 .67 k ( r) = single scattering albedo 


a = 0.87 exp (— 2k 4 r) + 0.97 [ 1 - exp (-2^)] = asymmetry factor 
kj = spectra! absorption coefficient of ice (Irvine and Pollack, 1968; Hobbs, 1974) 
r = radius of ice grains or the circular area equal to the shadow area of ice grains 
(Asano et al. 1979). 

Above equations constitute the model for the spectral incident flux under partially cloudy skies. 
The input parameters of the model are the clear sky atmospheric parameters (ozone path length, 
precipitable water, turbidity, and surface pressure), the cloud optical thickness and the grain size of 
snow. Illustrative results and comparison with observations follows. 


2-2. Illustrative Flux Results 

Since the emphasis in this report is or. the radiation characteristics over extensive snowfields, 
the atmospheric parameters representative of Antarctic coastal regions (such stations as Maudheim - 
71° 03' S 10° 56' W (Liljequist, 1956) and Mirny -66° 33' S 93" 01' l (Kusin, 1904)) was chosen: 
ozone path length 0.35 cm NTP (Robinson. I960). Angstiom turbidity 0.017 ( l.iliequist, 1^56), 


predpiteble water 0.25 cm (Nieman, 1977) and the surface pressure 1000 millibar. The ice grain size 
of snow is chosen to be 0.4 millimeters (Liljequist, 1956). The insolation (spectral integration from 
0.3 to 2.5 fa n) and the spectral flux values for a clear sky are shown in Figures 2*1 and 2-2 respectively. 
The observed range of insolation shown in Figure 2*1 is for Mirny during 1956 to 1957 period (Rusin, 
1964). The mean values of observed insolation at all solar elevations are same for Mirny and Maud- 
heim. The calculated insolation values differ from mean observations by about 10% at low elevations, 
ami about 6% at high elevations. These differences may be attributed to the inaccuracies in the 
parameterized transmission coefficients (Leckner, 1978) and in the estimation of atmospheric param- 
eters. The calculated values are definitely within the range of observations. 
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Figure 2-1 . Clear Sky Insolation 
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Figure 2-2. Clear Sky Spectral Flux at Different Solar Elevations 

For totally cloudy skies (c * 1 in Equation 2-1 ) the insolation values for different cloud optical 
thicknesses, together with observations at Mirny (Rusin, 1964) are shown in Figure 2 3. The mean 
values of the observed insolations are within 5 ' 7 < of the calculated values for the cloud optical thick- 
ness of 20. This optical thickness is observed to the mean for stratiform clouds < Feigelson. 1 ‘F’H I. 

The spectral flux values corresponding to this optical thickness is shown in Figure 2-4. A comparison 
of clear sky and totally cloudy sky spectral flux values. Figures 2-2 and 2-4, clearly illustrates the 
effect of multiple reflection between the cioudbusc and sr.owcover: the intensity of visible flux 
diminishes very little compared to n*'ar infrared llux. Polar ‘Wlnteouis' tCaichpok and Moodc . 

1971) is attributed to the intense and diffuse radiation which is observed <>• i : extensive snow Helds 
under cloudy skies. Table 2-1 gives the ratio of cloudy and clear sky total insolations for vatied clotiil 
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Figure 2-5: The Diffuse insolation As a Function of Fractional Cloud Cover 
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fraction and solar elevation. The daily average of this ratio for totally cloudy skies has been used in the 
parameterization of observed total insolation (Liljequist, 1956). The nonlinear variation of diffuse in- 
solation as a function of fractional cloud cover (Figure 2-5) shows clearly the effect of multiple reflec- 
tion (Schmetz and Raschke, 1 979). Comparison between calculated and observed insolations, shown 
above, illustrates that reasonable estimates of the input parameters of the model (atmosphere, cloud and 
snow) can provide good results for the fluxes. Computationally present model is little more expensive 
than the parameterized equations for the insolation (Sivkov, 1971; Lads and Hansen, 1 974) but it has 
the flexibility of providing the direct and diffuse spectral fluxes. 

The cloud model considered here is simplistic, partly because ante size of the cloud patches 
which leads to radiation leakage from the sides (McKee and Cox, 1974) is not considered. This re- 
finement, albeit important in any definitive study is difficult to incorporate into a simple model con- 
sidering the variability of cloud characteristics (shape and size of the cloud patches, and its distribu- 
tion). Present one parameter cloud model (the parameter being the optical thickness) is approximate 
but simple to be useful for practical applications e.g. in hydrology where an estimation of average 
effects of eloudeover may suffice. 

3. SPECTRAL REFLECTANCE OF SNOW 

The reflectivity of snow like many other natural media such as water and sand depends upon 
surface and subsurface radiative scattering and absorption. Snow is a highly complex and variable 
medium. Due to wind action and heat exchange with the atmosphere, the shape and size of ice crys- 
tals undergo considerable variation. If melt-freeze process occurs then irregular globular mass of 
polycrystalline grains as well as surface crust and internal ice lenses can form. A gale apart from 
compaction will renew the surface with fine grained rounded crystals. Accompanying these changes 
in the physical characteristics, are the changes in the reflectivity of snow. 

Observations with collimated radiation at m nlerate zenith angles (see Mellor, 1977) show that 
the re'lectance has a characteristic and almost p. ictable spectral dependence: it is high (generally 
greater than 907c) for visible radiation; with exception of resonance like structures near 1.1, 1.85 and 


2.2S pm the reflectance decreases to a few percent for near infrared radiation. In these observations, 
although the shape of reflectance curve did not change with time (except for the melt-freeze condi- 
tion) the magnitude at different wavelengths did change. With the aging of snow, the visible reflec- 
tance changed little compared to near infrared reflectance. 

For clear sky solar radiation, snow appeared to be a Lambertian reflector at high elevations but 
showed specular characteristics at low elevations (Knowles Middleton and Mungall, 1952; Dirmhim 
and Eaton, 1975). For overcaste skies snow behaved like a Lambertian reflector (Liljequist, 1956). 

Following previous studies on the reflectance of snow (Knowles Middleton and Mungall, 1952; 
Dunkle and Bevans, 1956; Liljequist, 1956; Bergen, 1975), a model for the spectral reflectance of 
snow was formulated in terms of specular surface reflection and volumetric multiple scattering 
(Choudhury, 1979c; Choudhury and Chang, 1980). For the surface reflection a bidirectional reflectance 
distribution function was defined which depended upon the crystal shape and topography. For volu- 
metric multiple scattering solutions of the radiative transfer equation were used which, for a deep 
homogeneous snowcover, depended upon the crystal size. For clear sky solar radiation, this model 
for spectral reflectance agreed with observed spectral and angular reflection characteristics. This 
comparison for albedo will also be shown in the next section. 

The spectral reflectance of snow can be calculated from the equation 

R(a) = p(a) f(a) + p(a) [ 1 - f(a)] Rj (a) + [ 1 - p(a)] R 2 (3-1 ) 


where 


p(a) - ratio of direct and total spectral fluxes for solar elevation a 
f(a) * directional surface reflectance 

w * 2 - 3q p 0 + (1 - «*) 3a*p 0 (2p 0 - q) 

Ri(a) , — 

2 (1 +q)(l -o 2 n 2 a ) 


(1 -a 2 ) w 
( 0 * * ... 

(1 - a 2 to) 


a* = a/(l + a) 





I 


o* (3(1 - w*) (1 - w*a*)] 1/2 

- 2p 
q “ 3(1 -w*a*) 

H 0 - sina 

and the remaining parameters are defined in Equation (2-7), which is used for the diffuse reflectance 
R 3 . The expression for the collimated reflectance R t , is obtained from the delta-Eddington approxi- 
mation (Joseph et al. 1976;Sobolev, 1975, Choudhury, 1979a;Wiscombe and Warren, 1979). 

The ratio of direct and total fluxes, p(a), is calculated from the incident flux model discussed 
in Section 2. The surface reflectance, Ua), calculated from the bidirectional reflectance distribution 
of a particulate medium (Trowbridge and Reitz, 1 975) is given in Table 3-1 for spherical, oblong 
and faceted surface crystals. Note that the surface reflectance of Bergen (1975) is defined through 
Fresnel reflectivity and specific surface, the latter being dependent upon density and grain size. 


Table 3-1 

The Directional Surface Reflectance for Different Grain Shape (Choudhury, 1979a) 



Solar Elevation (Degrees) 

I 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40.0 

Spherical 

0.337 

0.150 

0.089 

0.060 

0.044 

0.034 

0.027 

0.023 

Oblong 

0.360 

0.204 

0.139 

0.100 

0.073 

0.055 

0.043 

0.034 

Faceted 

0.660 

0.302 

0.194 

0.129 

0.086 

0.060 

0.044 

0.034 


*The grain shape is defined through an ellipsoidal representation. The ratio of major and minor 
axes for ‘descriptive* shapes spherical, oblong and faceted grains are, respectively, 1.0, 4.5 and 
31.5. Tire major axis is assumed parallel to the surface. The bidirectional reflectance distribu- 
tion function is from Trowbridge and Reitz (1 975). 

Few pertinent features of Equation 3-1 are worth noting. The spectral variation of reflectance 
arises primarily from the absorption coefficient of ice, which affects the single scattering albedo, u>, 
and the asymmetry function of scattering, a (Choudhury and Chang, 1979). Changes in the grain 
size affect the reflectance also through cc and a (see Equation 2-7). The contribution of surface re- 
flectance increases as the solar elevation decreases (see Table 3 1 ). At high elevations the spectral 
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reflectance and hence the albedo will vary primarily with the grain size. These dependencies will be 
dismissed Anther in the next section. 

Note that equation (3*1 ) is valid for a deep homogeneous snowcover. The effect of impurities 
and other internal inhomogeneities on the reflectance cannot be studied from this equation. Caution 
should be exercised when comparing the theoretical predictions with field observations. The visible 
reflectance is expected to be affected most by impurities and internal inhomogeneities due to large 
single scattering albedo and penetration depth. Changes in the surface grain, however, will affect the 
near infrared reflectance. 


4. ILLUSTRATIVE ALBEDO VALUES AND DISCUSSION 


The albedo of snow is defined as the ratio of reflected and incident solar energy. From Equa- 
tions (2*1 ) and <3-1 ) one obtains the albedo as: 


A(o,c) 


/R(a) G c (a) dX 
/G c (a)dX 


where the integration over wavelength (X) is from 0.3 to 2.5 pm, which encompasses most of the 
solar spectrum. The atmospheric parameters for the Antarctic coastal region, and typical stratus 
cloud optical thickness, discussed in Section 2, are used in the calculation of the incident flux. The 
ice grain radii chosen for the calculation of the spectral reflectance of snow are 0.1 5 and 0.4 milli- 
meters, which may be considered representative of spring and summer at the coastal regions (Lilje- 
quist, 1 956). The calculated albedo matrices (the albedo as a function of solar elevation and cloud 
fraction) are given in Tables 4*1 to 4-4. 

The albedo of snow generally shows diurnal variations. At Maudheim, the albedo of new fallen 
snow on December 2, 1951 was 0.77 which, during ‘warm’ weather, decreased to 0.73 on December 
8 but strong wind on December 1 1 gave 0.91 (Liijequist, 1956). It was argued that metamorphosis 
of the snow into larger rounded crystals during the day, and formation of crystals with more or less 
horizontal surfaces during the night occur. A gale, apart from compaction, renews the surface with 
fine grained and rounded crystals. At Mirny, the albedo decreased from 0.98 in July 1956 to 0.61 in 
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January 195 / (Rutin, 1964). A melt-freeze process leading to the formation of neve and ice was 
argued to be associated with this large change in albedo. Clearly these long-term changes in the albedo 
is associated with the metamorphism of snow i.e. changes in the crystal size and shape together wit!- 
surface texture. At Mirny, the mean solar elevation during July is much less than that during January, 
which also has contributed to the change in albedo, as discussed below. 

Apart from long-term diurnal variations the albedo depends upon solar elevation and cloudiness 
With respect to solar elevation, the albedo generally decreases as the elevation increases (Uljequist, 
1956; Rusin, 1964; Korff and Vonder H^r, 1972), but there is disagreement concerning the gradient 
of the albedo-elevation curve (see Korff «nd Vonder Haar, 1972, and Bryazgin and Koptev, 1969). 
Again due to metamorphic changes, the albedo is sometimes asymmetric with respect to the optimum 
elevation i.e. morning and afti moon values at the same elevation are not identical (Uljequist. 1956, 
Dirmhim and Eaton, 1975). The albedo under cloudy skies seems to depend upon the cloud altitude 
For high and mi ’die altitude clouds, ine albedo decreases as the solar elevation increases (Rusin, 1 9M > 
but for low altitude thick clouds the albedo tends to increase with the solar elevation ( Rusin, h»64; 
Korff and Von. 'or Haar, 1972). Compared to clear sky albedo, however, the albedo for cloudy skies 
is generally a tew percent higher. 

The calculated albedo matrices show tin. observed daily trend under clear skies and the effect 
of dease cloudct.ver For a clear sky, the v ;,r elevation dependence of the albedo depends upon the 
grain size. ..ml the shape of surface cij stai.,. At high elevations the albedo depends primarily upon the 
grain size. \t Sou elevation*. the albedo increases as the grains become faceted. The gradient of albedo 
elevation curve for sm. z v> : in; has undergone metamorphism under calm clear weatiier is expected to 
be higher than that for a wind blown snow. Associated with tins, one inuy also observe that at high 
elevations, the albedo of a wind blown snow is higher. Thus the magnitude and the gradient >>i .ubedo 
will depend upon the history of the snowcover. 

The spring (November) observations of albedo at Maudlteitn (l.iljcquist. 1 956 1 compare well with 
the calculated values (Figure 4-1 ). The mean gram size during spring is meaxiued to V ab. . 1 0 ? 
millimeters. The observed values are between those predicted tor rounded :::ul n -'at. |, 
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Figure 4*1 . Comparison of Observed and Calculated Gear Sky Albedo 


comparing the observations for December at Maudheim it was found that the grain size is about 0.4 
millimeters. Indeed, the measured grain sizes during spring and summer were respectively 0.2-0.4 mm. 
and 0.4-0.6 mm (Liifequtot, 1956V A comparison with observations at Mirny (Rusin, 1 964) is also 
shown in Figure 4-1 , The gradient of die albedo-elevation curve, both observed and calculated, is 
higher at Mirny compared to that at Maudheim. It is important to note that during observations at 
Maudheim red thawing was ohacrced only once, but significant metamorphic changes was observed 
at Mirny. Thus the physical characteristics of snow at these two stations were probably different 
during the observations. 

Barkstrom (1972), Mid Barks trom Mid Querfeld (1975) compared theoretical results with these 
observations without considering the spectral reflectnree and incident flux. Also these models do 
not consider the diffuse flux, which to of considerable importance at low elevations. Instead of 
calculating the albedo to terms of measured or expected grain size, an attempt war made to calcu- 
late it using other optical m ea sureme n ts (extinction coefficient and flux ratio). That being unsuccessful. 
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theoretical parameter! were limply adjusted to show “ . . . the computed albedos agree with ob- 
served albedos to within 1%, except at /u 0 ■ 0.08.” This approach to albedo is also taken by 
Wiscombe and Warren Q979). 

in this paper the abledo is not obtained by adjusting theoretical parameters, but as the ratio 
of reflected and incident energy using the spectral reflectance and the incident flux. l : or totally 
cloudy skies, the albedo tends to increase with solar elevation and it is a few percent higher than 
for clear sky values at least at high elevations. Quantitatively, the calculated values are in hettei 
agreement with observations at Maudhcim (‘about 90%,* Liljcquist 1956) than that at Mirny labour 
91%, Rusin 1964). It was shown by Bohren and Barkstrom (1974) that the spectral reflectance of 
mow for diffuse incident radiation is proportional to square root of the grain size. Present si mb i-i, ( 
the following approximate relation between the albedo and the grain size. 

A(a,c * I) *0.950 - 0.! 25 r 1 ” 

where r is the grain radius in millimeters. Because of inherent approximations, the formulation 
Bohren and Barkstrom (1974) cannot be used to calculate the albedo. 

A sense of dissatisfaction has been expressed regarding the ad hoc nature of separation of re- 
flected radiation into surface a«.J volumctrieally scattered components (Barkstrom, 1 9' 7 2). Ihe qu 
tion was posed as to whether the radiative transfer n snow should be treated as a problem similar m 
that of a plane parallel atmosphere or a water body since t )r the atmosphere, as opposed to a wat- 
body, this separation is not done. It is, however, difl'icm'. to identify snow with either of these tv, > 
media. Although for snow the surface is not a dielvctr c discontinuity like that of water or ice i ■■ 
difficult to visualize the ‘surface’ of atmosphere glittering like a snow' surface in clear sunny d ; i 
alizing that this separation is a debatable point, calculations were also done treat ino snow as a plain 
parallel atmosphere i.c f(a) B 0 in I quation 3-1 (figure 4-2). It can he seen that surface r. ! 1 ' i, iin . i 
the albedo at low elevations, and inclusion of this re lie*, tion provides letter agree m. r.t v.hli 

Sir. ,e an absolute accuracy oi the deltu-l-.ddmgtou approximation tor vodm, i tHf i 
tar has not yet been established for snow, we c<»mpatc.l the " subs of tins jppiovn.i"-y. ! i< 
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delta-two-stream (Schaller, 1979) and hybrid (Meador and Weaver, 19S0) approximations. The re- 
sults are shown in Figure 4-3. Difference between these approximations are quite noticeable for low 
solar elevations. Whereas the albedos for the delta-Eddington and the delta-two-stream fail to show 
observed inverse relationship with the solar elevation, the hybrid approximation does show this rela- 
tion. Note that inclusion of the surface reflection in the delta-Eddingto*; approximation essentially 
brings out this relation. Tire hybrid approximation, however, has been found to give inaccurate 
angular variation of albedo for thick peifect scattering media (see figure 4 in Meador and Weaver, 
1980). In the visible region snow is an almost pure scattering medium. It is therefore difficult to 
assess how truthfully the hybrid approximation is predicting the solar elevation dependence of snow 
albedo. The qualitative agree: v-nt between this approximation and observations of Liljequist (figure 
4-3) is to be noted, but without further study of this approximation it may not be used to draw con- 
clusions regarding the surface reflection. By parameter fitting Wiscombe and Warren (1979) errone- 
ously concluded that the delta-Eddington approximation without surface reflection can provide 
qualitative agreement with the observed solar elevation dependence of snow albedo. 

The model of snow considered here is simplistic. The shape and size of ice grains are c .-ife vari- 
able and the surface, being in direct contact with air, unc*; rgo metamorphism via wind action, vapor 
exchange and radiative processes. Observations led to stipulations that changes in surface texture and 
topography affect the albedo (Liljequist, 1956; Rusin, 1964; Hanson, 1960). Predictions from the 
radiative transfer equation will surely depend upon the shape and size of ice crystals (Cuzzi and 

Polalck, 1978), and the accuracy of the solution. This study is based upon approximate solutions of 
the radiative transfer equation. Obviously there is considerable room for further study. A ray tracing 
approach similar to Wolff (1975) should also be studied. For a better understanding of the reflective 
process, a careful set of spectral observations for wide elevation range is highly desirable. 
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Figure 4-2. Solar Elevation Dependence of Clear Sky Albedo 
for Different Crystal Shape 


5. CONCLUSION 

An approximate method for calculating the direct and diffuse spectral solar fluxes over extensive 
snowflelds for partially cloudy skies is discussed. The input parameters for the calculation of these 
fluxes are atmospheric precipitable water, turbidity, ozone content, surface pressure, the optical 
thickness of clouds and the grain size of snow crystals. A spectral reflectance model of snow was 
outlined which requires the grain size and grain shape as input parameters. Illustrative insolation and 
albedo values are obtained from spectral reflectance and incident flux for representative atmospheric 
and snow parameters of Antarctic coastal regions. Although there are uncertainties regarding the 
spectral values due to a lack of experimental results, the spectrally integrated values (insolation and 
albedo) compare favorably with observations. Unlike previous calculations of the albedo by ad- 
justing theoretical parameters, present study calculates it from the spectral values considering grain 





size and shape of ice crystals of snow. This study failed to reproduce the observed solar elevation 
dependence of tie albedo for dear sky conditions when radiative transfer in snow was treated as a 
problem similar to that of a plane parallel atmosphere. The conclusion and the agreement be- 
tween theory and observations are. however, tentative because the two-stream and the delta- 
aldington approximate solutions of the radiative transfer equation are used in this study. Fur- 
ther study of the radiative transfer equation for snow, and a a reful set of spectral reflectance 
measurement for a wide solar elevation range are highly desirable. 
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